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SUMMARY 


In  streamlined  flow  past  a  flat  plate  aligned  with  a  uniform  stream,  it  is 
shown  that  (a)  the  Goldstein  near-wake  and  (b)  the  Blasius  boundary  layer  are 
nonunique  solutions  locally  for  the  classical  boundary  layer  equations,  whereas 
(c)  the  Rott-Hakkinen  very-near-wake  appears  to  be  unique.  In  each  of  (a),  (b) 
an  alternative  solution  exists  which  has  reversed  flow  and  which  apparently 
cannot  be  discounted  on  immediate  grounds.  Thus,  depending  mainly  on  how  the 
alternatives  for  (a),  (b)  develop  downstream,  the  symmetric  flow  at  high  Reynolds 
numbers  could  have  2 ,  4  or  more,  simple  steady  forms. 

Concerning  non-streaml ined  flow,  e.g.  past  a  bluff  obstacle,  new  similarity 
forms  are  described  for  the  pressure-free  viscous  symmetric  closure  of  a  predom¬ 
inantly  slender  long  wake  beyond  a  large-scale  separation.  Features  arising 
include  nonuniqueness,  singularities  and  algebraic  behavior,  consistent  with  non¬ 
entraining  shear  layers  with  algebraic  decay.  Nonuniqueness  also  seems  possible 
in  reattachment  onto  a  solid  surface  and  for  nonsymmetric  or  pressure-controlled 
flows  including  the  wake  of  a  symmetric  cascade. 
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1 .0  INTRODUCTION 


The  first  part  of  this  work  is  concerned  with  a  class  of  basic  wake  flows 
which,  according  to  their  near-wake  properties  at  least,  are  found  to  possess 
nonunique  solutions  with  or  without  reversed  flow  present.  The  common  theme 
involved  is  the  nature  of  the  steady  laminar  planar  flow  past  a  finite  body, 
fixed  in  a  uniform  stream  for  instance,  at  large  Reynolds  number. 

At  one  extreme,  if  the  body  is  sufficiently  streamlined  then  attached  motion 
can  result  on  the  body  surface  (see  Section  5  however)  and  the  thin  wake  beyond 
is  governed  by  the  viscous  boundary  layer  equations.  At  the  other  extreme,  for 
non-streamlined  motion,  large-scale  separation  and  flow  reversal  take  place  in 
general  and  the  wake  is  thick  with  detached  viscous  shear  layers  present.  Even 
then  there  is  reason  to  think  that  the  main  recirculatory  eddies  may  well  be  so 
long  that,  on  the  characteristic  wake  scale  much  greater  than  the  body  dimen¬ 
sions,  the  wake  closure  and  subsequent  recovery  to  uniform  flow  occur  in  a 
predominantly  slender  flow,  i.e.,  the  wake  is  thick  but  very  long.  See  Sychev 
(1967),  Messiter  (1975),  Smith  (1979).  The  main  issue  there  centers  around  the 
flow  features  at  closure.  If  these  are  effectively  inviscid  then  vorticity 
conservation  tends  to  suggest  a  strong  reversed  jet  being  projected  upstream, 
along  the  center-line  in  symmetric  flow,  possibly  ruling  out  the  presence  of  only 
a  slowish  recirculatory  motion  overall.  So  the  alternative  of  significant 
viscous  action  being  brought  into  play,  before  or  after  the  main  wake  closure,  to 
suppress  a  strong  reversed  center-line  flow  upstream,  suggests  itself.  Thus,  for 
example,  the  main  closure  could  force  a  detached  shear  layer  to  return  almost  to 
the  center-line  by  inviscid  means  but  leave  any  necessary  reversed  flow  to  be 
sealed  off  by  the  viscous  wake  motion  holding  beyond.  The  above  two  extremes 
provide  the  motivation  for  the  study  of  viscous  wakes  under  certain  starting 
conditions.  Moreover,  in  both  extremes  the  induced  pressure  is  usually  negli¬ 
gible,  since  either  the  body  must  be  thin  to  be  streamlined  or  the  major  wake  is 
slender.  As  a  result,  although  there  are  some  exceptions  to  this  view,  one  of 
which  is  discussed  below,  we  are  led  to  consider  for  the  most  part  pressure-free 
viscous  wakes. 

Streamlined  motion  is  epitomised  by  the  symmetric  flow  past  a  finite  flat 
plate  at  zero  incidence.  The  pressure-free  boundary  layer  equations  control  the 
flow  on  and  beyond  the  plate,  apart  from  a  small  leading  edge  zone  and  a  triple¬ 
deck  at  the  trailing  edge  to  adjust  the  Blasius  solution  (see  Section  5)  to  the 
near-wake.  The  near-wake  here,  described  first  in  Goldstein's  remarkable  paper 
of  1930,  is  reconsidered  in  Section  2  below.  It  is  something  of  a  surprise  that 
the  solution  turns  out  to  be  nonunique:  one  solution  is  Goldstein's,  with  forward 
flow;  the  second  solution  has  reversed  center-line  flow.  Both  are  formally 
acceptable  as  far  as  we  can  tell,  with  the  second  solution  suggesting  an  account, 
for  the  trailing  edge  adjustment  and  the  downstream  wake,  quite  distinct  from  the 
well-known  classical  account. 
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Nonuniqueness  and  other  aspects  are  implied  also  for  nonsymmetric  or 
slightly  thicker  bodies  and  these  are  noted  in  Section  2.  In  contrast  the  Rott- 
Hakkinen  (1965)  "very  near"-wake  within  the  triple-deck  admits  only  one  solution 
for  symmetric  flow.  However,  the  extent  of  nonuniqueness  possible  in  the  aligned 
flat  plate  problem  is  enhanced  later  by  Section  5. 

The  application  of  pressure-free  viscous  wake  properties  to  non-streamlined 
flow  and  to  wake  closure  in  particular  is  discussed  in  Section  3,  by  means  of 
more  general  near-wake  forms.  Some  of  these  forms  lead  to  algebraic  singular¬ 
ities,  linking  with  the  author's  (1983)  proposal  of  shear  layers  exhibiting 
algebraic  decay  and  no  entrainment,  and  possibly  pointing  to  a  self-consistent 
account  of  large-scale  wake  closure.  Finally,  two  extra  kinds  of  near-wakes 
exhibiting  nonuniqueness  are  examined  in  Section  4.  One  concerns  the  long  wake 
behind  a  typical  body  in  a  symmetric  cascade,  where  again  the  boundary  layer 
equations  hold  but,  unlike  most  of  the  other  wakes  studied  here,  the  pressure  can 
be  significant  because  of  the  confinement  present.  For  a  wide  cascade  a  solution 
branching  from  the  uniform  flow  solution  is  found.  The  second  extra  near-wake 
considered  describes  the  behavior  of  slowish  eddy  motion  at  the  onset  of  wake 
closure. 

Secondly,  Section  5  concerns  flow  over  a  fixed  surface.  Aspects  similar  to 
those  of  Sections  3  and  4  apply  generally,  but  in  particular  it  is  conjectured 
that  even  the  pressure-free  boundary  layer  on  an  aligned  flat  plate  may  have  a 
nonunique  solution.  For,  as  well  as  the  Blasius  solution,  another  apparently 
self-consistent  starting  solution  near  the  leading  edge  is  possible  with  a 

reversed  flow  coming  to  rest  there.  In  all  these  examples  of  new  solutions,  * 

however,  it  remains  to  be  seen  whether  some  external  agency  must  be  present, 

e.g.,  downstream,  to  allow  a  deviation  from  the  usual  solutions,  or  not.  If  so, 

the  nonuniqueness  found  could  still  be  of  physical  significance.  But  if  no 

external  agency  is  required,  then  in  view  of  the  possible  nonuniqueness  of  the 

near-wake  (Section  2)  as  well  as  the  Blasius  solution  (Section  5)  there  could  be 

2,  or  even  more,  simple  steady  solutions  for  that  most  basic  problem  of  free  ♦ 

symmetric  flow  past  a  flat  plate  in  a  uniform  stream. 

Throughout  these  near-wake  and  boundary  layer  flows  the  planar  boundary 
layer  equations  apply,  usually  pressure-free  and  hence  classical,  although 
similar  features  including  nonuniqueness  can  be  expected  to  persist  in  nonzero 
pressure  fields  and  in  three-dimensional  motions.  Further  details,  especially  of 
the  contexts  from  which  this  work  arose,  including  non-entraining  shear  layers 
and  the  cascade  flow  problem,  are  given  by  Smith  (1982a,  1983).  Other  contexts 
include  the  motion  of  liquid  layers  and  the  flow  near  or  above  a  finite  heated 
plate,  for  example.  The  nondimens ional  scaled  variables  u,  J>,  x,  Y,  p  are  used 
below,  denoting  the  streamwise  velocity,  the  stream  function,  the  streamwise  and 
lateral  coordinates,  and  the  pressure,  respectively,  and  Re  is  the  Reynolds 
number.  We  have  in  mind  subsonic/incompressible  fluid  flow,  although  the 
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implications  for  supersonic  motion  are  similar.  All  the  cases  of  nonuniqueness 
and/or  reversed  flow  and/or  singularities  found  below  in  near-wakes  and  on  solid 
surfaces  pose  quite  challenging  coraputat ional  tasks,  classical  or  interactive, 
for  the  rest  of  the  wake  or  boundary  layer  or  for  the  shorter  scale  flows 
implied.  It  would  seem  worthwhile  addressing  these  tasks,  nevertheless,  e.g. 
to  test  the  nonuniqueness  of  flow  past  a  flat  plate,  the  extent  of  that  non¬ 
uniqueness,  and  to  continue  the  study  of  large-scale  wake  closure.  These 
problems  remain  of  both  fundamental  concern  and  wide  application. 


*Also,  with  nonuniqueness  present,  the  stability  of  the  flows  and  their 
unsteady  development  from  an  initial  state  assume  extra  significance. 
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2.0  THIN  NEAR-WAKES 


The  thin  viscous  wake  of  a  flat  plate  (y  ■  0,  -l<x<0)  placed  at  zero 
incidence  to  a  uniform  stream  is  governed  by  the  classical  boundary  layer 
equations  (with  y  *  Re-1/2Y) 

u«li  >ulu  -  M  (2.1a,b) 

3Y  3x  3x  3Y  3Y* 

holding  in  x>0,  with  the  boundary  conditions 


u(x,«)  -  1,  ij>(x,0)  *  —  (x»0)  -  0,*  (0+,Y)  -  ♦b(Y).  (2.  lc-f) 


Here  ’f'g(Y)  in  (2. If)  is  Che  Blasius  profile  (see  also  Section  5)  at  the 
trailing  edge  x  *  0“,  satisfying  tyB(0)  “  ♦«'(())  “  0  (no  slip). 

*B"(0)  -  A(»  0.33206...),  4»b* (-)  -  1  and  2*B'"  +  <f»B*B”  *  °*  while  (2-lc-e) 
are  the  free  stream  and  symmetry  conditions. 

The  near-wake,  for  x  -*■  0+,  then  has  the  Goldstein  (1930)  double  structure, 
to  accommodate  the  nonlinear  change  from  no  slip  to  symmetry  at  Y  ■  0,  Across 
most  of  the  near-wake  the  expansion 

+  -  *B(Y)  +  x1/3  Aj  tl»B’(Y)  +  ...  (2.2a) 


holds,  for  0<Y<«,  describing  a  displacement  of  the  oncoming  Blasius  flow 
controlled  by  the  unknown  constant  Aj .  Near  Y  ■  0,  in  contrast,  \|i  takes  the  form 

♦  -  x2/3  G(£)  +  ....  with  l  =  Y/x1/3  ,  (2.2b) 


because  of  (2.1e),  and  from  (2.1a,b,d-f)  G  satisfies 
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—  GG”  -lG'2  -  0 


3  3 

G(o)  =  G"(o)  =0 


G(5)  ~  i  X  \2  +  o  (C2)  as  I 
2 


(2.3a) 

(2.3b,c) 

(2.3d) 


The  constraint  (2.3d)  affords  the  match  with  (2.2a)  as  well  as  (2.1f),  since 
(2.3a)  allows  G'~X£  +  c  with  the  constant  c  unknown,  so  that  X  A,  ■  c. 

The  similarity  problem  (2.3a-d)  has  a  solution  due  to  Goldstein,  with 


G'(o)  =  1.611  X2/3,  giving  c  =  0.892  X2/3  ,  (2.4) 


which  then  fixes  Aj  ■  X-1c. 

The  main  issue  of  this  section,  however,  is  that  the  solution  is  not  unique. 
A  second  solution  of  the  nonlinear  problem  (2.3a-d)  exists  which  may  be  equally 
acceptable  on  present  grounds,  i.e.,_in  the  context  of  (2.1a-f).  The  second 
solution  arises  as  follows.  Taking  G‘(o)  to  be  nonzero,  since  otherwise  G(£)  =  o 
from  (2.3a-c),  we  set 


G  -  |G* (o) | 1/2  G(£),  l  -  |G'(o)|~1/2  C.  (2.5) 

which  casts  the  problem  for  G(£)  as  an  initial  value  one: 


G'"  ♦  I  G  G"  -  i  G'2  -  o,G(o)  -  G”(o)  -  o,G’(o)  -  ±  1,  (2.6a-d) 

3  3 


where  ±  correspond  to  G'( o)  ^  0,  from  (2.5).  Also,  the  condition  (2.3d)  now 
becomes  an  unknown  asymptotic  property. 


G~—aU  +  b)2, 
2 


(2.7) 
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to  be  acquired  by  the  solution  of  (2.6a-d),  it  is  hoped,  with  G'(o),  c  then  being 
found  from  |G'(o)|  *  (X/a)2/3  and  c  =  a1/3  b  X2/3.  See  also  Weyl  (1941,  1942a, b) 
Jones  &  Watson  (1963). 

Analysis  as  in  Section  3  below,  supported  by  a  numerical  treatment  of 
(2.6a-d)  (Fig.  1),  shows  that  both  values  *1  for  G'(o)  yield  the  behavior  (2.7). 
If  G'(o)  ■  1  we  obtain 


a  **  0.489  ,  b  =  1.132  ,  (2.8a) 

giving  the  original  Goldstein  results  (2.4);  but  if  G'(o)  =  -1  we  find 

a  *  6.082  ,  b  »  -3.048  ,  (2.8b) 

which  give  the  alternative  results 

G'(o)  »=  -0.300  X2/3  ,  c  =  -5.564  X2/3  (2.9) 


for  the  solution  of  (2.3a-d).  This  second  solution  has  the  properties  that 
G'(Ci)  *  0,  G(C2)  ■  0  where  £,  •  2.129,  £2  4  3.376.  Hence  reversed  flow  (u<o)  is 
present  for  0  Y  ]g*(o)[— 1/2  x1/3,  near  the  wake  center-line,  as  x  +  0+ , 

and  the  separating  streamline  1(1  ■  o  is  at  Y  *  (2  |G’(o)|-1/2  x1/3.  The  center- 
line  velocity  is  reversed,  u(x,o)  ■  -0.300  X2/3  x1/3  for  small  positive  x. 

The  two  solutions  are  shown  in  Figs,  l(a-d).  Figure  1(a)  shows  these  solu¬ 
tions  in  terms  of  the  normalized  velocity  functions  of  Eqs.  (2.6a-d)  for  the 
region  near  the  wake  centerline.  The  same  results  are  also  shown  in  Fig.  1(b)  in 
terms  of  the  velocity  functions  of  Eqs.  (2.2b),  (2.4)  and  (2.9).  Also  shown  in 
Fig.  1(b)  are  the  rapid  approach  to  the  asymptotic  solution,  Eq.  (2.3d).  No  more 
(regular)  solutions  can  exist,  incidentally,  from  the  above  argument.  Neither  do 
solutions  exist  for  any  negative  oncoming  shear  X. 

It  is  intriguing  to  consider  the  above  nonuniqueness  in  relation  to  the  wake 
flow  further  downstream  and  to  the  flow  nearer  the  trailing  edge.  If  the  first 
solution  (2.4)  of  Goldstein  holds  then  the  wake  motion  remains  entirely  forward 
for  all  x  >  0  and  can  be  determined,  necessarily  numerically,  by  forward  marching 
in  x  from  the  starting  forms  (2.2a,b)  with  (2.4).  Far  downstream,  as  x  ♦  •  ,  the 
uniform  state  is  achieved  in  the  form 
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♦  ~  Y 


n 


Y/x 


1/2 


(2.10a) 


consistent  with  the  momentum  deficit  integral 


00 

I0  “  J  u(l-u)dY  =  constant  (>  o) 
o 


(2.10b) 


holding  throughout  x  >  o,  from  (2.1a-f),  and  determined  by  the  Blasius  form  of 
(2. If).  The  motion  is  also  forward  in  the  triple-deck  and  other  regions  induced 
closer  to  the  trailing  edge  where  x  is  0(Re~3/8)  and  smaller;  see  Jobe  &  Burggraf 
(1974)  and  references  in  Messiter  (1979),  Stewartson  (1981),  Smith  (1982b).  An 
apparently  complete  description  of  the  flow  past  the  aligned  flat  plate  at  large 
Re  is  achieved,  and  comparisons  with  experiments  and  numerical  Navier-Stokes 
solutions  are  favorable  at  surprisingly  low  Reynolds  numbers. 

If  the  second  solution  (2.9)  holds,  on  the  other  hand,  then  reversed  flow 
u  <  o  must  be  present  for  at  least  an  0(1)  distance  x  in  the  wake  as  depicted  in 
Fig.  1(c).  The  solving  of  (2.1a-f),  again  necessarily  numerical,  is  a  more 
difficult  task  due  to  the  upstream  influence  in  the  reversed  motion.  The 
integral  (2.10b)  is  preserved  still,  suggesting  (2.10a)  far  downstream,  but  here 
and  below  the  question  of  regularity  at  wake  closure,  where  u(x,0)  0  with 

x  -*■  xQ  >  0,  also  has  to  be  addressed  (indeed,  could  there  be  more  than  one 
closure  point,  implying  multiple  eddies  in  the  wake?).  Local  analysis  indicates 
that  a  close  analogue  of  Goldstein's  (1948)  singularity  for  instance  cannot  occur 
at  a  wake  closure  here,  especially  with  negligible  pressure  gradient,  but  a 
regular  description  is  not  likely  either.  For  (2.1b)  and  its  Y-derivatives  then 
suggest  that  only  an  identically  zero  profile  is  allowed  at  x  =  Xq.  With  a 
pressure  gradient  acting  regular  or  singular  wake  closure  can  take  place,  e.g., 
analogous  to  Goldstein's  (1948)  forms,  A  prime  example  of  a  wake  with  non- 
uniform  pressure  occurs  in  a  branching  flow  where  the  flat  plate  upstream  is 
sited  within  a  channel,  parallel  to  the  walls:  Bates  (1978),  Badr  et  al  (1983). 
Without  a  pressure  gradient  however  it  may  well  be  that  one  of  the  singularities 
in  Section  3  below  is  encountered  as  x  ♦  xQ  -,  or  xQ  may  be  infinite.  If 
there  is  a  non-removable  singularity  at  closure  that  would  tend  to  suggest  that 
an  external  agency  is  required,  downstream  perhaps,  to  allow  (2.9)  to  occur  in 
the  near-wake,  and  the  free  flow  past  a  flat  plate  adheres  to  the  original  form 
(2.4).  If  not,  however,  (2.9)  introduces  nonuniqueness.  Removal  of  any  singu¬ 
larity  then  presumably  would  invoke  an  interactive  local  account  near  x  ■  x0. 

The  second  solution  (2.9)  must  also  be  able  to  serve  as  a  downstream  (X  ♦  «) 
condition  for  the  triple-deck,  the  most  important  of  the  flow  structures  nearer 
the  trailing  edge,  in  which  x  *  Re“3/8X.  There  the  central  point  is  that  the 
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problem  of  adjusting  from  the  Blasius  flow  (as  X  ♦  -*»)  to  the  near-wake  (as 
X  ■+  +»)  is  elliptic  in  character,  because  of  the  local  pressure-displacement 
interaction.  Hence  there  appears  to  be  no  obvious  reason  to  discount  the  second 
solution  (2.9)  as  a  downstream  asymptote  there  (Fig.  1(b)):  it  is  consistent 
with  the  controlling  equations  at  least.  On  both  scores,  then,  for  smaller  or 
larger  x,  the  second  solution  (2.9)  does  not  seem  inconsistent  yet.  It  would  be 
interesting  to  conduct  computational  studies  of  the  two  new  nonlinear  problems 
posed  by  (2.9):  the  triple-deck  and  the  wake,  both  having  reversed  flow  of 
course.  Supposing  that  each  problem  does  have  a  solution  then  we  have  apparently 
a  second  description  of  the  flow  past  an  aligned  flat  plate  at  large  Re.  If  so, 
the  nonuniqueness  would  be  expected  to  persist  at  finite  Re  and  Navier-Stokes 
calculations  then  would  be  of  interest.  Again,  nonuniqueness  could  continue  to 
hold  as,  say,  slight  degrees  of  nonsymmetry  are  introduced  into  the  flow  field. 

An  alternative,  however,  is  that  the  triple-deck  problem  admits  the  solution 
(2.9)  as  a  downstream  state  only  when  the  local  trailing  edge  motion  is  suffi¬ 
ciently  disturbed,  say  by  a  small  nonzero  incidence  or  body  thickness  as  in 
Melnik  &  Chow  (1975),  Ruban  &  Sychev  (1979),  Smith  &  Merkin  (1981),  Smith 
(1982c).  As  the  disturbance  increases  in  size,  could  the  second  solution  (2.9) 
(with  reversed  flow)  come  into  play  before  or  after  separation  occurs  at  the 
trailing  edge,  for  instance?  Could  (2.9)  enable  the  range  of  triple-deck  solu¬ 
tions  with  separation  to  be  enlarged?  In  all  such  cases  it  is  necessary,  for  the 
above  nonuniqueness,  that  the  wake  problem  (2.1a-f)  holding  in  x  >  0  has  a  solu¬ 
tion  with  the  starting  form  (2.9)  and  that  should  be  regarded  as  a  primary 
concern.  Section  5  below  is  also  relevant  here,  for  the  flow  ahead  of  the 
trailing  edge. 

Another  feature  concerns  the  "very-near-"  wake  within  the  triple-deck 
itself,  as  X  •*  0+.  There  a  local  similarity  form  holds  again  but  governed  by  the 
Rott-Hakkinen  (1965)  problem 

I'”  8  8"  -  j  I'2  -  |  Pi  (2.11a) 

g(o)  *  g"(o)  -  0,  i'  -  XT  l  *  0  as  l  *  •  (2.11b-d) 


in  symmetric  flow,  with  the  constant  >  0.  Here  the  local  wake  has  to  be 
displacement-free,  from  (2. lid),  and  the  pressure  coefficient  pj  is  unknown,  in 
contrast  with  (2.3a-d).  Ihe  solution  of  (2.11a-d)  is  found  to  be  unique  (with 
g'(o),  pj  taking  Hakkinen  &  Rott’s  (1965)  values),  however,  according  to  a 
normalization  based  on  g'(o)  as  before  and  integration  for  p1  varying.  Figure 
1(d)  gives  the  displacement  function 
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A  =  “lira 
a  -  E,  -*•  oo 


£  -  - 


I"(") 


as  it  varies  with  p^  for  g'(o)  ■  ±1  as  well  as  the  asymptotic  solution  for 
pj  ♦  «°,  corresponding  in  effect  to  g’(o)  ■  0.  Note  that  only  one  solution  gives 
A  *  0  as  required  by  Eq.  (2. lid)  and  that  is  Rott  and  Hakkinen's  (1965).  Also, 
no  solution  producing  g''(®)  finite  exists  to  the  left  of  the  dotted  curves, 
which  both  give  |a|  ♦  •  at  finite  negative  values  of  pj .  Thus  nonuniqueness  is 
not  compounded  by  the  very-near-wake ,  in  symmetric  flow  at  least. 
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3.0  THICKER  NEAR- WAKES 


Thicker  bodies  or  larger  angles  of  incidence,  i.e.  non-streamlined  motions, 
tend  to  produce  large-scale  separation  and  eddies  eventually  (e.g.,  Cheng  &  Smith 
1982)  and  there  the  properties  of  wake  closure  downstream  have  to  be  addressed. 

We  consider  symmetric  flow,  for  which  the  governing  equations  near  the  center- 
line,  beyond  the  major  wake  closure,  are  again  the  classical  boundary  layer  ones. 
If  in  addition  the  main  eddies  upstream  are  predominantly  slender,  as  in  the 
extended  Kirchhoff  approach  for  example,  then  the  induced  pressure  is  small, 
being  proportional  to  the  eddy  slope,  and  so  the  pressure- free  equations  (2.1a,b) 
apply.  Also,  from  an  origin  shift,  x  >  0  again,  and  (2.1c-e)  hold.  Only  the 
starting  profile  u(0+,Y)  at  x  *  0+  remains  to  be  specified.  If  it  has  forward 
velocities  only,  smooth,  nonzero  and  of  order  unity,  then  (2.1a-e)  can  be  inte¬ 
grated  forward  in  x,  giving  a  fairly  standard  wake  phenomenon.  If,  on  the  other 
hand,  the  wake  closure  is  not  quite  complete  yet,  then  the  starting  profile 
contains  some  reversed  velocities  and  the  wake  is  non-standard.  We  consider 
first  whether  (a)  such  reversed  flow  can  be  present  near  the  center-line  and  (b) 
it  can  be  initially  small.  Here  (a)  is  meant  to  correspond  to  the  downstream 
condition  for  a  smaller  scaled,  mainly  inviscid,  process  of  near-closure  in  which 
for  instance  an  incoming  detached,  forward  moving,  shear  layer  is  joined  to  a 
reversed  center-line  motion.  Condition  (b)  then  allows  that  reversed  motion  to 
remain  weak  upstream  and  hence  consistent  with  a  slowish  (e.g.,  extended 
Kirchhoff)  eddy  there.  So  the  idea  is  that  complete  closure  is  achieved  through 
the  action  of  viscosity  downstream  when  x  is  finite,  whereas  a  complete  closure 
for  small  x  would  be  inviscid  and  force  a  strong  reversed  center-line  motion 
inconsistent  with  the  slowish  eddy  ahead.  Further  details  here  are  given  by 
Smith  (1983).  Non-standard  wakes  can  arise  also  if  (c)  the  starting  profile 
gives  forward  but  small  flow  anywhere. 

The  above  and  other  contexts  indicate  a  local  viscous  similar 
holding  in  the  near- wake  as  x  0+  for  some  Y  >  0  of  the  form  \|>  “ 
with  C  ■  Y/x*^K+*  .  Here  K  is  a  constant  with  |k|  >  1  so  that  u 
G'(£)  is  small:  see  later  for  | K |  <  1.  From  (2.1a,b,d,e)  the  function  G 
satisfies 


ity  solution 

xR/<K+1>GU) 

,  X(K-1 ) / (K+l  ) 


(K+l  )G' ' '  +  KGC"  -  (K-l)G'2  -  0  (3.1a) 


G(o)  -  G' '(o)  ■  0. 


(3.1b,c) 
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Also,  from  the  reasoning  in  Section  2  (the  special  case  K  *  2)  and  a  normaliza¬ 
tion,  two  extra  conditions  are  allowable, 


G'(o)  -  ±  1. 


(3. Id) 


Hence  the  possibility  of  nonuniqueness  arises  for  any  starting  profile  of  the 
kind  (a)-(b),  or  (c). 

Computational  solutions  of  (3.1a-d)  are  presented  in  Figs.  2(a-e).  They 
fall  into  line  with  certain  analytical  properties  below  and  with  the  view  that  as 
£  increases  from  zero  only  two  ultimate  trends  can  be  observed.  If  K  <  -1  a 

singularity  occurs  at  a  finite  position  £  *  £0  >  0,  with 

C  ~  C !  ( C0  -  C)“  1K  I  +  0(£o  -  O'1,  £  +  £Q  -.  (3.2a) 


In  contrast,  if  K  >  1  the  solution  continues  for  all  £,  giving  (apart  from  an 
origin  shift) 


G  ~  C2  £K  +  0(£-1) ,  £■*•«.  (3.2b) 

Here  Cp  C2  are  unknown  constants  but  both  are  positive,  for  an  integral  of 
(3.1a-c)  gives 

G"  -  /« «,)«,}/  =,2<«2>“p{jcr/2<:<tl>d«l}«2 

(3.3a) 


which  shows  that  the  vorticity  is  positive. 


G"  >  0  ,  for  £  >  0, 


(3.3b) 


since  |k|  >  1.  Therefore  the  velocity  G'(£)  increases  with  £,  from  the  center- 
line  values  in  (3. Id),  after  which  the  ultimate  trends  (3.2a,b)  both  produce 
forward  flow,  G'  ♦  ♦■».  It  is  interesting,  as  seen  in  Fig.  2(a),  that  these 
trends  are  inevitable  whether  G'(o)  is  positive  or  negative,  for  that  yields 
nonuniqueness  as  the  trends  (3.2a,b)  represent  features  of  the  given  starting 
profile. 


■  V 

!> 
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Both  (3.2a,b)  form  self  consistent  asymptotes  for  (3.1a),  but  if  K  <  -1  the 
asymptote  G  ~  C2  £K  as  £  -*■  »  is  unobtainable  as  it  requires  a  change  in  sign  of 
the  vorticity,  contradicting  (3.3b).  A  contradiction  also  arises  then  f.-om 
substitution  into  (3.3a).  Similarly  the  asymptote  G  ~  Cj(£0-£)K  as  (  +  £0- 
is  unobtainable  for  K  >  1.  Two  other  asymptotes  acceptable  for  (3.1a)  alone  are 
G  ♦  constant  G(«)  as  (  +  »  and  G  ~  6(£-£0)_1  either  as  (  +  £0-  or  as  (  +  «. 

In  the  first  however  G(«°)  must  be  positive,  to  avoid  exponential  growth  as  £  -*■  «, 
and  that  contradicts  (3.3b).  In  the  second,  the  vorticity  is  negative  if  £0  is 
finite,  while  it  must  change  sign  if  £  ♦  «,  so  that  again  neither  limit  is 
acceptable  for  the  whole  problem  (3.1a-d).  Substitution  into  (3.3a)  again 
confirms  this.  This  second  asymptote  corresponds  to  an  exact  solution  G  * 
6(£-£0)-1  of  the  governing  equation  (3.1a),  incidentally:  see  Smith  (1983), 
Appendices  A,  B,  and  Section  4(ii)  below. 

For  K  >  1,  then,  (3.2b)  is  achieved.  Also  the  viscous  region  here  (£  “ 

0(1))  is  small  and  contracting  as  x  ♦  0+ .  Outside,  for  Y  of  0(1),  the  starting 
profile  has  the  form  ^(O^.Y)  ~  C2Y^  as  Y  +  0+,  to  match  (3.2b),  and  as  shown  in 
Fig.  2(b,c)  overall  the  flow  structure  and  nonuniqueness  are  similar  to  those  of 
Section  2.  Some  limiting  analysis  is  possible:  Appendix  C. 

The  properties  implied  for  K  <  -1,  where  the  viscous  region  is  large  and 
expanding  as  x  ♦  0+,  are  quite  distinct.  There  the  singularity  (3.2a)  occurs 
(see  also  Appendices).  So  near  £  ■  £„  a  thinner  region  arises,  far  from  the 
center-line,  wherein  Y  -  £Qx-*^  |K|~iy=  ?  is  Qf  0(1)  and_the  starting  profile 
<»(0+,Y)  is  arbitrary  except  that  4»(0+,Y)  ~  CjjY|  as  Y  +  -  •»,  to  match  (3.2a), 
and  a  matching  condition  is  required  as  Y  •*•  «.  The  latter  depends  on  how 
multi-structured  the  near-wake  is:  see  Figs.  2(d)  and  (e).  The  simplest  condi¬ 
tion  is  3i|)(0+,  ®)/3Y  »  1,  which  satisfies  (2.1c)  directly  and  gives,  for  the 
wake  displacement  6(x)  defined  by  i|i  ~  Y  -  6(x)  as  Y  ♦  the  growth 


6(x)  ~  £q  x-1/<lKl"1>(>  «),  x  ♦  0+. 


(3.4) 


This  growth  upstream  seems  sensible  physically  with  regard  to  both  the  wake 
closure  anticipated  further  downstream  and  to  the  detached  flow  properties 
holding  on  shorter  scales.  Also  the  algebraic  decay  (?♦-<■)  of  the  0(1) 
thick  region  of  faster  flow  above  is  in  line  with  the  algebraic  decay  attainable 
in  detached  shear  layers  on  a  shorter  length  scale:  see  Smith  (1983).  Again, 
other  realistic  conditions  are  possible  for  Y  ♦  •.  A  decay 
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*(0  +  ,  b  ~  C3  ♦  C4  Y~\  i  **  -,  (3.5a) 

like  that  for  Y  +  -  ®,  can  be  tolerated  with  v(>  0),  C4 ,  C3  (a  mass  flux) 
constants.  That  induces  a  further,  thicker,  viscous  region  outside  of  thickness 

Y  ■  0(x~*/^~^)  and  tji  -C3  ■  )  is  small,  provided  v  >  1.  In  effect 

(3.1a)  holds  again  with  K  replaced  by  -v  and  (3.1b-d)  are  replaced  by 


G  ~  C4U  “  ?0)"v  as  C  ♦  S0  ♦  , 


(3.5b) 


from  (3.5a).  As  £  increases  further,  then  slower  motion  results, 


G  ~  6£-1  or  constant,  as  £  *  <*>, 


(3.5c) 


consistent  with  (3.1a);  or  a  singularity  like  (3.2a)  is  repeated,  and  so  on.  See 
also  Appendices.  An  analogue  of  (3.3a)  can  also  be  applied.  Conversely,  if 
V  <  1  then  no  such  outer  region  arises  and  a  direct  join  to  a  further  incoming 
shear  layer  more  removed  from  the  center-line  can  be  made.  Thus  the  solutions 
with  K  <  -1  seems  perhaps  more  relevant  to  a  shorter  scale  detached  flow  than 
those  with  K  >  1,  if  significant  reversed  flow  is  to  proceed  into  the  viscous 
wake . 


Other  points  are  as  follows.  First,  simple  positive  power  solutions  for 
G(£)  in  (3.1a)  exist  only  for  K  *  1,2;  see  also  (4.7a,b)  below.  For  K  m  1, 

G(£)  *  ±  constant  £  satisfies  (3.1a-c)  but  gives  u  *  ±  constant,  producing  a 
center-line  flow  not  coming  to  rest,  in  contrast  with  a  slowish  eddy  motion  in 
the  shorter  scale  features.  For  K  ■  2,  G(£)  ■  -dj£2  satisfies  (3.1a,b)  but  not 
(3.1c).  Instead  a  no  slip  condition  is  satisfied,  a  matter  taken  up  in  Section 
5.  Second,  a  range  of  solutions  of  (3.1a),  but  with  no  slip  G  ■  G'  ■  0  at  (  “  0, 
exists:  again  see  Section  5  below.  Third,  if  a  significant  pressure  gradient, 
classical  or  interactive,  is  acting  then  more  possibilities  open  up:  e.g.,  there 
are  Kennedy's  (1964)  classical  solutions;  Gajjar  &  Smith  (1983)  demonstrate 
branching  in  liquid  layer  and  hypersonic  flows;  and  ultimate  forward  flow  as  in 
(3.2a,b)  can  be  avoided  in  special  cases,  e.g.,  Goldstein's  (1948).  The 
relevance  to  external  wake  closure  behind  a  slender  eddy  is  then  apparently  lost, 
however.  Fourth,  alternative  accounts  of  the  near-wake  are  obtainable  of  course 
with  the  center-line  motion  reversed  but  strong,  i.e.  u(0+,  Y)  finite  and 
arbitrary  for  finite  Y,  with  0(0*,  0)  negative  say.  The  eddy  flow  ahead  must 
then  contend  with  a  strong  reversed  flow,  possibly  by  means  of  algebraic  decay 
again  and  non-entrainment  (Smith  1983).  Fifth,  6(»)  »  I0  from  (2.10b).  Sixth, 
an  integral  form 
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3u 

3Y 


(3.6) 


of  (2.1a-e)  with  V  *  -3<|»/3x  confirms  the  property  (3.3b)  that  3u/3Y  >  0  for  small 
Y  >  0  if  the  center-line  flow  is  coming  to  rest,  U(x,o)  ♦  0,  as  x  ♦  0'f.  Hence 
the  velocity  u  has  a  local  minimum  at  the  center-line.  Seventh,  as  in  Section  2, 
numerical  studies  of  the  complete  wakes  stemming  from  the  above  starting  forms 
with  or  without  reversed  flow  would  be  of  interest,  as  would  continued  studies  of 
the  shorter  scale  detached  motions.  Eighth,  the  energy  integral  of  (2.1a-e) 
requires 


—  1  5  —  I  /  u(l-u2)dY} 
3x  1  E)  3x  l  ’ 


(3.7) 


So  due  to  dissipation  the  energy  thickness  6£(x)  increases  monotonical ly  with  x 
after  usually  being  positive  at  the  start  of  the  wake.  Finally,  the  similarity 
solutions  above,  and  below,  apply  also  to  flow  at  x  ■  0“ ,  coming  to  or  from  rest 
if  |k|  >  1.  There  |x|  replaces  x,  so  that  the  nonlinear  terms  in  (3.1a)  change 
sign.  Then  letting  G  ♦  -G  reproduces  (3.1a-d),  hence  (3.2a,b)  and  hence  the 
conclusion  that  sufficiently  far  from  the  center- line  the  motion  must  be 
reversed,  whether  the  center-line  velocity  is  positive  or  negative.  This  con¬ 
clusion  can  be  avoided  in  special  cases,  e.g.,  Goldstein's  (1948),  only  if  a  non 
xero  pressure  gradient  is  present. 


♦ 
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4.0  ADDITIONAL  NEAR-WAKES  AND  NONUNIQUENESS 


Two  other  kinds  of  near-wakes  of  interest,  particularly  concerning  wake 
closure  and  the  possibility  of  nonuniqueness,  arise  in:  (i)  the  wake  of  a  non- 
streamlined  body  placed  in  a  wide-spread  cascade  or  wind-tunnel;  (ii)  the  slowish 
flow  between  a  detached  shear  layer  and  a  wake  center-line.  Their  properties  are 
summarized  below. 

(i)  The  wide-spread  cascade.  Consider  the  boundary  layer  equations 

„  -  !i  ,  u  2H  -  iH  -  -  i£  +  (4 . la) 

3Y  3x  3x  3Y  dx  3YZ 


governing  the  0(Re)  long  wake  of  a  body  forming  part  of  an  infinite  cascade 
(Fig.  3),  such  that  the  symmetry,  mass  flux  and  center-line  conditions 


^  »  in 

3Y 


0  at  Y 


o;  ♦  -  l.  in 

3Y 


0  at  Y 


(4.1b,c) 


hold.  Here  i^,  x,  Y  have  been  normalized  with  respect  to  the  0(1)  upstream  mass 
flux,  to  Re,  and  to  the  0(1)  cascade  spread,  and  unlike  in  most  of  the  other 
wakes  studies  in  this  paper  the  pressure  p  is  unknown  here.  For  a  relatively 
wide-spread  cascade  or  thin  body  the  starting  condition  is  simply  the  original 
uniform  stream 


u  *  1,  i|>  -  Y  at  x  »  0+,  for  0  <  Y  <  1  (4.  Id) 

c.f.  a  f inite-spread  cascade  would  have  u  ■  (1-H)~*  for  H<Y<  1,  u  *  0  for 
0  <  Y  <  H,  at  x  *  0+,  with  H  between  0  and  1. 

One  solution  of  the  wake  problem  (4.1a-d)  is  the  undisturbed  original  flow 
u  -  1,  *  ■  Y,  p  «  0,  for  all  x  >  0,  Y,  (4.2) 


of  course,  meaning  that  wake  closure  and  recovery  have  taken  place  already  before 
the  final  stage  of  (4.1a-d)  is  entered.  There  is  another  possibility,  however, 
suggested  by  an  alternative  near-wake  form  of  (4.1a-d),  as  follows. 
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Suppose  the  flow  near  the  center-line  is  small,  locally  uniform  and  reversed 
near  x  *■  0,  given  by  -u  =  0(xn),  *  0(Yxn),  n  unknown.  To  adjust  the 

incoming  stream  (4. Id)  to  this,  for  small  x,Y,  a  thin  shear  layer  is  necessary, 
of  the  Chapman  form  and  therefore  requiring  an  entrainment  *  0(x1/2).  Hence 
the  reversed  motion  below  the  shear  layer  can  supply  the  entrainment  if  the  shear 
layer  is  positioned  at  Y  *  0(xl/2~n),  requiring  n  <  1/2.  The  majority  of  the 
near-wake  above  the  shear  layer  then  must  accommodate  the  0(xl/2-n)  displace¬ 
ment  but  preserve  mass  flux,  from  (4.1c),  and  so,  effectively  from  one-dimen¬ 
sional  channel  flow  theory,  a  pressure  drop  0(x1/2-n)  is  induced.  This  agrees 
with  the  0(x2n)  pressure  fall  necessary  to  drive  the  reversed  center-line 
motion  if  n  *  1/6.  Formally  therefore  we  have  the  expressions  for  small  x, 


i|>  ~  -YjY  x*/6  ,  u  ~  -Yj  x1^  ,  for  0  <  Y  <  y2  x^3  (4.3a) 


below  the  shear  layer,  where  Ylt  Y2  are  unknown  positive  constants;  the  corre¬ 
sponding  pressure 


P 


-1/2 


Y?  xl/3  ; 


(4.3b) 


the  Chapman  shear  layer 


*  x1^2fc(n)  +  ...,  u  -  f^.(n)  +  ....  for  Y  -  Y2x1/^3  ♦  *1/,2n  »  (4.3c) 

where  2f£"  +  fcf^.'  “  °»  *£(")  “  1.  *£(“•)  “  °»  8ivin8  “  “< 

with  tc  -  1.24;  and  finally 

Y  2  2 

♦  *  Y  +  x1/3  — (Y-l)  +  ....  u  -  1  +  x1^3  — i—  +  ...  (4.3d) 

2  2 


above  the  shear  layer,  to  conserve  mass.  All  of  (4.3a-d)  satisfy  the  governing 
equations  (4.1a-d)  for  small  x.  Matching  (4.3a,c)  and  then  (4.3c,d)  gives 
YLY2  -  K,  Y2  -  Yi2/2  in  turn,  yielding 
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Yj  “  (2<) 1 ,  Y2  “  <  <2/2 ) 1 


(4.3e) 


The  local  expansion  then  appears  to  be  self-consistent  and  further  terms 
including  exponentials  can  be  generated  in  principle  (SmitK  ’983).  The  starting 
form  (4.3a-e)  represents  a  second  ( eigen) solut ion,  branching  irom  the  incoming 
uniform  stream  and  an  alternative  to  the  first  solution  (4.2).  Note  that  the 
momentum  integral  equation 


p(x)  + 


/  u2dY  -  1 
o 


(4.4) 


obtained  from  (4.1a-d)  is  satisfied  by  (4.3a-e).  Indeed,  there  appears  to  be  no 
theoretical  objection  to  the  second  account  (4.3a-e)  at  first  sight,  nor  any 
obvious  reason  why  regular  wake  closure  and  subsequent  flow  recovery  back  to  the 
uniform  state  (4.2)  should  not  take  place  further  downstream,  for  x  >  0,  yielding 
self-consistency  as  seen  in  Fig.  3.  Here  the  branching  solution  for  cascade  flow 
is  depicted  showing  the  small-x  layers  0(x1/3),  0(x1/2)  flow  reversal  and  the 
possible  eddy  closure  downstream.  Accordingly  it  would  be  very  interesting  to 
see  the  results  of  a  numerical  study  of  the  branching  solution,  with  reversed 
flow,  stemming  from  (4.3a-e).  The  alternative  to  the  structure  depicted  in 
Fig.  3  is  simply  the  original  uniform  stream,  u*l ,  for  all  x. 


(ii)  Flow  between  a  detached  shear  layer  and  a  wake  center-line.  If  such  flow 
is  sufficiently  slow  then  it  too  is  controlled  by  the  viscous  boundary  layer 
equations  but,  unlike  most  other  motions  considered  in  this  paper,  it  can  also  be 
affected  by  the  small  pressure  gradient  due  to  a  predominantly  slender  eddy.  For 
the  flow  to  be  so  slow  the  entrainment,  into  the  shear  layer  which  shelters  the 
slow  flow  from  faster  motion  outside,  must  be  obliterated  since  otherwise  the 
slow  moving  fluid  is  effectively  inviscid.  The  mechanism  for  non-entrainment 
involves  algebraic  decay  in  the  shear  layer  (Smith  1983)  and  hence  an  algebraic 
growth  within  the  slower  flow,  wherein  we  have  in  scaled  terms 


u 


3u  _  j3u 
3x  3x  3Y 


— p 1  ( x  ) 


(4.5a,b) 


*  «  .  0  at  Y  -  0  , 

3Y 


(4 ,5c ,d) 
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if  there  is  no  center-line  jet.  Here  S(x)  gives  the  shear  layer  position,  and 
for  0  <  Y  <  S(x).  The  algebraic  growth  requires 


6(x-xj) 
(Y-S(x) )  * 


6(x-xj  ) 
(Y-S(x))^ 


as  Y  +  S(x)- , 


(4 . 5e) 


consistent  with  (4.5a,b)  and  with  the  shear  layer  (Smith  1983).  The  constant  Xj 
is  unknown. 

We  seek  local  solutions  of  (4.5a-e)  for  eddy  closure,  as  x  ♦  o  ±  say, 
supposing  that  the  slowish  flow  here  can  speed  up  then,  i.e.,  |u]  -*•  ®.  Thus 
*  ~  xK/(R+l)G(g)t  y  -  6  x1/(K+1)  again,  u  ~  X(K-1)/(K+1)  c.(e)  but  now  |K|  <  lf 

for  x  •*  0+.  So  G  satisfies  (3.1a-d)  if  the  pressure  gradient  becomes  negligible 
then.  Numerical  solutions  for  |k|  <  1  are  shown  in  Fig.  4,  while  (3.3a)  now 
requires 


G"  <  0  for  5  >  0. 


(4.6) 


As  a  result  the  ultimate  trends  encountered  are  not  covered  by  (3.2a,b).  Instead 


G  ~  6U  -  ?0)_1  as  K  *  $0  -  ,  (4.7 a) 

or  G  ♦  G(»)  >  0  as  ?  ♦  -  ,  (4.7b) 

or  G  ~  C2  +  0(T1)  a*  5  ♦  •  .  (4.7c) 


Here  (4.7a)  is  attained,  for  some  £0  >  0,  whenever  G'(o)  *  -1;  (4.7a)  also 
holds  if  G'(o)  ■  1  when  -1  <  K  <  but  for  <  K  <  K2,  G'(o)  «  1,  (4.7b)  is 
attained,  while  (4.7c)  is  found  for  G' (o)  ■  1 ,  K2  <  K  <  1 .  Here  Kj,  K2  ■  0.5  and 
are  almost  certainly  equal,  a  result  suggested  by  analysis  in  Appendix  B. 

Further  analytical  support  can  be  given,  e.g.,  Appendices  A-C.  Mote  that  the 
limits  G  ~  65"*,  5  +  •,  and  (3.2a)  are  not  attainable. 


The  algebraic  growth  anticipated.  (4.5e),  (4.7a),  can  be  achieved  for  all 
|k|  <  1,  then,  giving  S(x)  ~  50x*''R+1*  tending  to  Tero  as  at  a  closure  point  at 
x  ■  0+ .  Closure  at  x  ■  0“  has  |x|  replacing  x  above,  G  changes  sign  in  (3.1a) 
and  again  (4.7a)  can  be  achieved. 


Alternative  local  solutions  are  possible  with  a  significant  pressure 
gradient  (Appendix  C),  with  ♦  ,  u  regular,  or  with  Xj  zero  (or  x  ♦  x^l).  Further 
details  concerning  (i),  (ii)  above  are  given  by  Smith  (1983). 
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5.0  NONUNIQUENESS  OF  THE  BLASIUS  SOLUTION, 
AND  SURFACE  REATTACHMENT 


This  second  part  of  the  report  comes  as  a  by-product  of  Section  3,  stemming 
from  the  exact  solution  G  *  -djC2  of  (3.ia)  with  reversed  flow,  for  K  *  2,  which 
gives 


u  *  -2d jY ,  =  -djY2  . 


(5. la) 


The  constant  d^  is  assumed  positive  here.  The  local  solution  (5.1a)  for  small  x 
satisfies  the  following  conditions.  (a)  No  slip  holds  at  Y  ■  0.  (b)  The  Chapman 

entrainment  i|>  =  -xx1/2  is  supplied  at  the  position 


Y  -  (.e/dj)172  x  1/4  +  .... 


(5.1b) 


around  which  the  thinner  0(x1/2)  Chapman  form  (see  Section  4(i))  applies.  (c) 

The  reversed  velocities  u  induced  between  Y  *  0  (the  solid  surface)  and  the  shear 
layer  at  (5.1b)  are  small,  of  order  x1/4  as  x  ♦  0,  in  keeping  with  a  slowish  eddy 
flow  upstream,  (d)  The  pressure  response  required  remains  negligible. 

The  small-x  form  (5.1a,b)  as  depicted  in  Fig.  5  satisfies  the  boundary  layer 
equations  (2  la-c)  with  (2.1d-f)  replaced  by 


^=ua!0  at  Y  e  0, 


(5.2a) 


u  “  1  at  x  ■  0  for  Y  >  0  . 


(5.2b) 


Further  terms  in  the  small-x  expansion  stemming  from  (5.1a,b)  (see  also  Appendix 
D)  can  be  generated  at  will,  integral  properties  can  be  derived,  and  once  again 
it  would  be  interesting  to  follow  the  subsequent  downstream  development,  of 
(2. la-c)  with  (5.2a,b),  which  necessarily  is  a  numerical  and  challenging  task 
since  reversed  flow  is  present. 
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The  starting  form  (5.1a,b)  and  its  induced  displacement 


6(x)  ~  ( </d i )  1  / 2  x1/4  ,  x  ♦  0+ 


(5.3) 


describe,  first,  what  is  in  effect  (from  2.1a-c),  (5.2a,b))  the  flow  over  an 
aligned  flat  plate  in  a  uniform  stream,  and  so  provide  an  alternative  to  the 
Blasius  solution  ~  xi/2  ♦jj'fh),  n  “  Y/x*/2  for  which 


6(x)  *  x1/2  ,  x  +  0+  . 


(5.4) 


Hence  the  suggestion  is  that  even  this  most  basic  flow  problem  may  exhibit  non¬ 
uniqueness.  If  so,  then  for  some  d1  (5.1a)  -  (5.3)  must  constitute  also  a  multi- 
structured  downstream  asymptote  for  a  smaller  0(Re“2/®)  interactive  zone  close  to 
the  leading  edge,  again  an  alternative  to  the  Blasius  asymptote  for  the  more 
orthodox  Navier-Stokes  0(Re-1)  zone.  Whether  (5.1a)  -  (5.3)  can  apply  for  the 
whole  flat  plate,  or  whether  an  external  agency  is  required  further  downstream 
say,  remains  unknown  as  yet,  and  the  question  of  a  removable  or  nonremovable 
singularity  downstream  also  arises,  despite  a  Goldstein  (1948)-like  one  being 
ruled  out  by  the  comments  following  (5.5a-c)  below  and  as  in  Section  2.  However, 
the  above  suggests  a  novel  form  of  leading-edge  separation  and,  coupled  with  the 
suggestions  of  Section  2,  there  is  now  the  chance  of  2,4  or  even  more,  simple 
solutions  being  found  for  the  high  Reynolds  number  steady  flow  past  a  finite  flat 
plate  at  zero  incidence. 

The  second  main  application  of  (5.1a)  -  (5.3)  implied  is  to  the  start  of 
viscous  pressure-free  reattachment  onto  a  solid  surface,  downstream  of  the  gross 
separation  and  eddy(ies)  produced  by  an  obstacle  on  the  surface  or  by  a  ramp  for 
instance.  The  starting  profile  can  be  more  general  than  (5.2b),  incidentally, 
provided  it  is  positive  at  Y  *  0+.  As  in  the  neaT-vakes  before,  a  pressure- 
displacement  interaction  does  not  have  to  be  present,  neither  is  a  strong  center- 
line  flow  produced  (from  (5.1a,b)).  The  induced  pressure  affects  the  flow 
significantly  only  within  a  shorter  length  scale  and  then  only  near  the  surface 
at  first  where  the  reversed  eddy  motion  is  slowest,  from  (5.1a,b).  There  (Smith 
1983)  the  possibility  of  secondary  separation  occurs,  depending  on  the  interac¬ 
tion  present. 


When  compared  with  Section  2  there  is  an  origin  shift  x  ♦  x  -  1. 
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The  simple  starting  form  (5.1a)  is  in  fact  a  special  case  of  the  solutions 
of  (3.1a)  but  with  no  slip  holding, 


G(o)  =  G'(o)  =  0  (5.5a,b) 

G"(o)  =  ±  1.  (5.5c) 


Comments  similar  to  those  near  (3.1a-d)  apply  here,  with  ip  =  XK/(K+1)G(£)  and 
Y  *  X1/(K+1)?.  Solutions  are  given  in  Fig.  6  for  Eq.  (3.1a)  with  the  constraints 
of  (5.5a-c),  corresponding  to  eddy  closure  on  a  solid  surface,  as  K  varies  with 
G''(o)  =  ±1.  For  our  main  interest,  i.e.  fluid  coming  to  rest  as  x  +  0*,  |k|  >  1 
again  and  then  (3.2a,b)  are  repeated  when  K  <-l ,K  >  1  in  turn,  with  but  one  range 
of  exceptions:  K  =  2  for  which  (5.1a)  holds  (G  =  ±1/2  £2),  allowing  forward  or 
reversed  flow,  and  the  cases  1  <  K  <  2  with  G''(o)  *  -1  which  give  (4.7a)  and 
ultimate  reversed  motion.  All  other  cases  by  contrast  give  forward  flow 
ultimately,  from  (3.2a,b),  including  K  =  3;  c.f.,  a  nonzero  pressure  gradient 
would  allow  Goldstein's  (1930),  (1948)  separation/reattachment  singularity  to 
occur  for  K  =  3  with  reversed  flow  if  x  +  0+  or  forward  flow  if  x  +  0"  (see  Fig. 
6).  So  a  generalization  of  the  Chapman  match  (5.1b)  seems  unlikely,  although 
further  properties  akin  to  those  of  Section  3  are  found.  See  also  Appendices  A- 
C. 


As  in  Section  4(ii),  the  solutions  with  |k|  <_  1  have  a  different  application 
concerning  eddy  closure  and  different  features.  Solutions  are  shown  in  Fig.  6, 
comments  similar  to  those  of  Section  4(ii)  apply  and  some  properties  can  be 
obtained  from  Appendices  A-C. 

As  in  Section  3,  the  corresponding  similarity  properties  of  flow  over  a 
solid  surface  as  x  +  0-  can  be  deduced  by  change  of  sign  of  G,  yielding  ultimate 
reversed  flow  if  |k|  >  1  except  for  the  range  1  <  K  <  2  described  above. 
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APPENDIX  A 

SPECIAL  CASES  OF  NEAR-WAKES 


Helpful  results  for  the  near-wake  similarity  form  (3.1a-d)  are  obtained  in 
certain  special  cases. 

First  we  consider  K  ♦  1  ±  with  G'(o)  -  -*-1 .  Then  G  *  q£  +  eG1  +  0(e2)  if 
K  x  1  +  e.  Here  q  =  1,  |e|  is  small  and  Gj  satisfies  a  linear  equation  with 
Gj(o)  »  G^o)  *=  Gj''(o)  ■  0  from  (3.1a-d).  The  solution  has 


Gj"  .Iq2e-qC2M  J5eqCi2/4d£i  #  (A1 ) 

^  o 

so  that  Gj ' 1  ~  £-1  as  £  *  ®.  Hence  Gj  ~  £  In  £,  giving  G  ~  £1+E,  £  +  in 

line  with  (3.2b)  for  K  *  1+  and  (4.7c)  for  K  =  1“. 

Second,  suppose  K  ♦  1±  but  with  G'(o)  =  -1 ,  a  less  simple  case.  For  £  of 

order  unity  G  expands  as  above  but  with  q  =  -1 ,  so  that  (A1 )  gives  exponential 

growth, 


"  ~i  wl/2e52M  as  ^  +  . 


(A2) 


A  nonlinear  phase  enters  for  large  £,  therefore,  with  £  «*  a  +  a“*£  where  the 
equation  |e[  exp  (a2/A)  =  a3  determines  positive  a  »  1,  %  is  0(1)  and  G  =  a  8 
with  6  of  0(1).  Hence  S(?)  satisfies  the  Chapman/Bias ius  equation,  from  (3.1a), 


28"  '  ♦  8  8"  -  0, 


(A3) 


with  6 - 1  ±  2x^2e^2  as  t  ♦  -«*,  (A4) 

for  e  <  0,  to  match  (A2).  Here  (A3),  (A4)  pose  a  shooting  problem  for  8.  If 
e  >  0  the  perturbation  in  (A4)  is  positive  and  (A3),  (A4)  yield  in  effect  the 
Chapman  solution  for  which  8' (-•)  ■  0  and 
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2' +  constant  2'(*°)  >  0  as  %  *  «  . 


Hence  2  ~  a22'(“)  (5  -  a)  for  £  >  a  to  leading  order,  in  keeping  with  (3.2b) 
again.  If  c  <  0  the  perturbation  in  (A4)  is  negative  however  and  then,  since 
2' '  <0,  (A5)  is  unobtainable.  Instead  2  acquires  the  counterpart  of  (4.7a)  at  a 
finite  value  of  %• 

Third,  if  K  =  0  then  (3.1a-d)  reduces  to  solving 


H"  +  H2  =  0 


H(o)  *  ±1,  H’(o)  =  0 


(A7,A8) 


for  H  =  G' .  Here  (A6)  gives  3H' 2  +  2H3  =  ±  2  on  integration,  so  that  since 
H*  <  0 


-(ip-  i 


d  h 

(±l-h3)1/2 


Therefore  (4.7a)  is  obtained  (H  ♦  -<*>)  in  both  cases,  with 


1/2C  =  f  _ d_h__ 

\  3/  °  i,  (il-h3)1' 


(A10) 


fixing  the  values  of  £0  in  turn.  This  agrees  with  Section  4(ii)  and  Fig.  4. 

A  Fourth,  suppose  K  ■  -1  +  e  with  |ej  small.  Then  the  viscous  region  thins, 

K  ■  being  the  characteristic  scale,  and  G  ~  |c|^2G.  So  we  obtain  from 

(3.1 a-d ) 


±G' ' '  -  G  G"+  2  G|Z  ■  0 


(All) 


G(o)  -  G"  (o)  -  o,  G'(o)  -  ±  1, 


(A12) 
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leading  50  the  singularity  G  *  | K0  ~  C|  1ln  |  £Q  -  as  ?  +  £  -  for  some  positive 
CQor  as  ?  ♦  in  every  case.  The  thinning  and  the  singular  response  here 
reflect  the  behavior  of  the  solutions  of  Sections  3  and  4(ii)  for  K  near  -1.  In 
addition  the  case  (All),  (A12)  corresponds  to  an  expfmjx)  (|x|  -*•  <*)  or  exp(-m2/ 

| x ] m3 )  (jx|  *  0)  form  of  x-dependence  for  <i,  Y  with  m^,  m2,  m3  being  constants. 

Fifth,  for  large  |k]  the  ultimate  form  is  initially  exponential,  in  line 
with  the  asymptotic  injection  problem  (c.f.  Jones  &  Watson  1963),  before 
algebraic  behavior  as  in  (3.2a,b)  takes  over. 

Perturbations  about  the  two  solutions  of  Section  2  for  K  ♦  2  ±  also  confirm 
the  trends  of  (3.2a,b)  holding  them. 

If,  on  the  other  hand,  a  pressure  gradient  is  present  then  (3.1a-d)  becomes 
more  like  the  Falkner-Skan  problems.  Some  solutions  with  bounded  G'(»)  exist  as 
in  Jones  &  Watson  (1963),  Kennedy  (1964),  while  Gajjar  &  Smith  (1983)  note  solu¬ 
tions  with  algebraic  growth.  Special  cases  can  also  be  analyzed  as  in  (A1 )  - 
(A12). 
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APPENDIX  B 

OTHER  ASPECTS  OF  THE  SIMILARITY  PROBLEMS 


There  is  some  value  in  regarding  the  similarity  problems  (3.1a-d)  and  others 
in  a  phase  plane,  with  G  versus  G' .  For  the  possible  forms  of  ultimate  behavior 
of  G  from  the  differential  equation  (3.1a)  itself  can  be  listed  as  (3.2a)  (K  <  -1 
only),  (3.2b)  (K  >  -1  only),  (4.7a),  (4.7b)  and 

G  ~  6  £-1  as  £  +  «  ,  (Bl) 


and  it  is  useful  to  decide  which  forms  can  be  achieved  from  the  boundary  condi¬ 
tions  (3.1 b— d ) . 

If  | K |  >  1  as  in  Section  3  then  (3.3b),  coupled  with  G'=  dG/d£,  implies  that 
as  £  increases  the  solution  of  (3.1a-d)  must  move  to  the  right  in  the  G'-  G 
plane,  and  upward  in  the  right  half  (G1  >  0)  and  downward  in  the  left  (G1  <  0)  as 
depicted  in  Fig.  Bl .  Here  the  typical  phase  plane  G'-  G  for  equations  (3.1a-d) 
is  shown  when  | K ]  >  1.  The  asymptotes  are  designated  by  solid  arrows  la,  lb,  2, 

3,  and  4  while  the  open  arrows  denote  the  general  direction  of  the  solutions  in 

view  of  equation  (3.3b).  Hence  for  K  >  1  and  G'(o)  *  ±  1  the  forms  (3.2a), 
(4.7a,b),  (Bl)  are  unattainable  and  instead  (3.2b)  is  inevitable.  This  includes 
the  case  K  *  2  of  Section  2.  For  K  <  -1  however  only  (3.2a)  can  be  reached, 
given  the  movement  above. 

If  | K |  <  1,  on  the  other  hand,  as  in  Section  4(ii),  (4.6)  applies  instead 
and  the  movement  for  increasing  £  must  be  to  the  left,  although  upward /downward 
as  before.  So  all  of  (4.7a,b),  (3.2b)  are  attainable  in  principle  for  G'(o)  =  1, 
but  only  (4.7a)  is  accessible  for  G'(o)  ■  -1 ,  as  in  Fig.  B2 .  Here  the  phase 
plane  G'-  G  is  depicted  for  equations  (3.1a-d)  when  |k|  <  1  and  G'(o)  ■  ±  1. 

Again  the  open  arrows  give  the  general  directions  for  increasing  £  and  the  solid 
arrows  denote  the  limits.  Note  that  infinitely  many  solutions  approach  the 
limits  la  and  3.  For  G'(o)  “  1  limiting  analysis  (Appendix  A)  shows  that  (4.7a) 

holds  for  K  *  0  but  (4.7c),  i.e.,  (3.2b),  for  K  ■  1-.  Between  them  there  is  a 

separatrix,  believed  to  be  for  the  value  of  K  in  (B4)  below,  at  which  (4.7b) 
holds  since  above  this  value  (4.7a)  is  inaccessible.  Perturbations  about  (B4) 
below  behave  similarly  to  (A3)  -  (A5)  above  and  support  the  response  near  the 
separatrix. 

A  similar  view  in  a  phase  plane  can  be  adapted  to  verify  the  properties 
achieved  in  Section  5  and  Fig.  6,  for  the  no  slip  conditions. 
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Another  aspect  concerns  the  limiting  forms  (Bl),  (4.7a),  the  exact  solution 
G  *  6£_1  of  (3.1a),  and  their  accessibility.  With  the  origin  shift  in  5  absorbed 
if  necessary,  we  set  G  -  6C*1  small  and  then  (3.1a)  shows  that  generally 


G  ~  65-1  +  aU|q  +  o(U|q)  , 


(B2) 


for  |£j  ♦  0  or  |C|  ♦  “,  where  (K+l)q3  +  3(K-l)q2  +  2(4K-5)q  +  12  K  *  0.  There¬ 
fore  either  q  *  -2 ,  corresponding  again  to  an  origin  shift,  or  (K+l)q2  +  (K-5)q  + 
6K  *  0,  giving 

q  -  (5  -  K  ±  (25-34K  -  23K2 ) 1 /2| /  (2R  ♦  2)  .  (B3) 


There  are  usually  two  roots  for  q.  They  are  both  real  if  25  -  34K  -  23K2  >  0, 
i.e.,  for  K“  <  K  <  K+  where 


K1 


-17 


23 


0.53886  ... 
-2.01725  ... 


(B4  ) 
(B5) 


but  the  roots  are  complex  conjugates  if  K  <  K-  or  K  >  K+  in  which  case 


real  (q)  -  (5  -  K)/(2K  +  2) 


(B6) 


from  (B3).  The  roots  or  their  real  parts  are  given  in  Fig.  B3.  Because  of  (B2) 
the  £  -  large  asymptote  (Bl)  is  possible  only  if  q  <  -1  or  real  (q)  <  -1 ,  in 
(B3),  (B6)  respectively,  and  this  restricts  K  to  the  range  -7  <  K  <  -1 .  Outside 
that  range,  for  K  <  -7  or  R  >  -1 ,  the  C  -finite  asymptote  (4.7a)  is  a  candidate 
instead.  This  holds  apart  from  the  special  values  K  *  K*  where  (B3)  has  equal 
roots,  suggesting  inaccessibility  in  line  with  Fig.  4.  Also,  for  K  ♦  -1  the 
logarithmic  behavior  noted  in  Appendix  A  takes  over. 

We  note  that  provided  -7  <  K  <  -1  it  is  possible  for  a  solution  of  (3.1a)  to 
proceed  smoothly  from  the  form  (3.2a),  at  (  ■  £e-  say,  to  the  form  (B2)  for 
C  ♦  -“,  as  C  decreases,  according  to  (B3),  (B6)  and  the  phase  plane.  This 
property  is  used  by  Smith  (1983)  in  the  nonuniqueness  of  shear  layers. 
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APPENDIX  C 

SPECIAL  CASES  OF  NO-SLIP  FLOWS 


These  arise  in  the  similarity  forms  (3.1a)  with  (5.5a-c)  considered  in 
Section  5.  C.f.  Appendix  A. 

First,  suppose  K  -*•  2  ±  with  G''(o)  -  +1.  Letting  K  -  2  +  e  with  |e|  small, 
we  have  G  *  1/2 qC ^  +  eGj  +  0(e2)  where  q  =  1.  Hence  (3.1a),  (5.5a-c)  yield  the 
solution. 


with  px  »  0  here.  So  G1"'  ~  £-1  as  (  +  »,  giving  Gj  ~  1/2  £2  In  £  and  G  ~  1/2 
£2+e,  in  line  with  (3.2b)  for  K  *  2 ±. 

Second,  if  K  =  2  +  e  with  |e|  <  <  1  again  but  G''(o)  *  -1,  then  q  =  -1  above 
and  (Cl)  gives 

Gj'"  ~  I  exp  (£3/9)  as  £  ■►  »  (C2) 


where  I  ■  3-2/3  (—1/3)1  >  0.  The  growth  (C2)  leads  to  a  new  nonlinear  stage  at 
large  £  *  b  +  b-2  X,  with  ?  of  0(1),  |e|  exp  (b3/9)  *  b8  fixes  b  »  1, 
and  G  "  b2  £.  There  the  shear  layer  problem 


'  +  1 
3 

5  5"  -i  g'2  -  0  , 

3 

(C3) 

5  ~  -  1/2  ± 

27  I  exp  (?/3)  as  t  ♦  , 

(C4) 

holds,  from  (3.1a),  (C2),  for  e  ^  0.  If  e  >  0  the  positive  perturbation  in  (C4) 
is  consistent  with  the  shear  layer  solution  in  which  &'(-")  «  0,  5' '(+«)  >  0  (see 
Stevartson  1981,  Smith  1982),  yielding  G  ~  (l/2)b6C' '(•»)( £-b)2  for  £  >  b,  in 
agreement  with  (3.2b)  for  K  ♦  2*.  But  if  e  <  0  the  negative  disturbance  in  (C4) 
causes  S  to  become  singular  at  a  finite  value  of  ?  in  agreement  with 
(4.7a).  Both  these  cases  agree  with  the  properties  of  Fig.  6  near  K  ■  2  when 
G"(o)  ■  -1. 
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Third,  we  consider  K  *  0,  for  which  (3.1a),  (5.5a-c)  require  (A6)  but  with 


H(o)  -  0  ,  H'(o)  ■!!,  (C5.C6) 

where  H  =  G' .  Hence  3H' 2  +  2H3  *  3,  and  if  H'(o)  *  -1  we  have  H  <  0  and 


H 

C  -  -  f  - 4*2 -  ,  (C7) 

o  (l-2/3h3)1/2 


giving  the  singularity  (A. 7a)  with  H  ♦  and 


£  «  f  _ ™ _ 

-«  ( l-2/3h  3) 1/2 


(C8) 


If  H'(o)  ■  +1 ,  however,  H  ^  0  for  0  £  £  £  2£j  and  there  is  symmetry  about  £  *  £j. 
In  0  <  £  ^  (j  the  first  sign  in  (C7)  in  reversed  and  H'  0.  Here 


(3/2) 


1/3 


/ 


dh 


(1-2/3  h3)1/2 


At  £  ■  £j,  H  attains  its  maximum  (3/2)1/3,  while  in  £  £j 


H 

f  dh 

(3/2)1/3  ( l-2/3h  3) 1/2 


(C9 ) 


(CIO) 


So  H  passes  through  zero  at  £  ■  2£j  before  acquiring  the  singular  behavior  (A. 7a) 
again,  with 


*o  "  <1  + 


(3/2) 1/3 

/ 

— «D 


dh 

( l-2/3h3) 1/2 


(Cl  1 ) 


The  results  (C8),  (Cll)  are  in  line  with  those  of  Fig.  6  for  K  •  0. 
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Comments  concerning  K  -*•  -1  and  j K |  ♦  •»,  as  well  as  pressure  gradients,  are 
similar  to  those  of  Appendix  A.  Note  that  as  in  Gajjar  &  Smith  (1983)  the  pres¬ 
sure  term  p^  in  (Cl)  for  instance  can  be  adjusted  to  avoid  the  growth  (C2)  and 
hence  avoid  (4.7a)  if  necessary.  A  like  effect  can  be  achieved  in  (Al),  (A2)  and 
elsewhere.  Finally,  perturbations  about  the  two  solutions  for  K  ■  1,  one  the 
Blasius  and  the  other  yielding  (4.7a),  confirm  the  occurrence  of  (4.7a,c)  for 
K  -*>  1±. 
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APPENDIX  D 

FURTHER  TERMS  IN  THE  ALTERNATIVE  FLAT  PLATE  SOLUTION 


Of  the  further  terms  generated  in  the  small-x  expansion  (5.1a,b)  the  most 
"dangerous"  and  interesting  are  exponentially  small  ones  within  the  reversed 
motion.  Algebraically  small  ones  overall  pose  little  difficulty.  The  exponen¬ 
tials  stem  from  the  behavior  of  the  Chapman  solution  at  its  lower  extremes, 


fc(n)  ~  -<  +  Ae<r^2 


(Dl) 


as  n  =  [Y-(x/dj )  l/2x1/l*]x_1/2  ♦  -®,  where  A  is  a  positive  constant.  As  a 
result  the  slower  reversed  motion  between  the  Chapman  layer  and  the  surface  must 
have  the  underlying  development 


♦  ■  -d^x1/2  +  ...  +  xMA(x)  exp  [g(x)/*1^]  .  (D2 ) 


Here  M  is  a  constant,  A,g  are  unknown  functions  of  x>  but  it  is  required  that  the 
Real  (g)  <  0  for  exponential  decay  as  x  ♦  0+,  and  x  =  x“1/l*Y  with  0  <  x  < 
(</dx)1/2. 

Substituting  (D2)  into  the  governing  equations  (2.1a,b)  we  obtain  the 
differential  equation 


x(g  ♦  Xg')  *  g'2 


(D3) 


for  g(x)»  where  without  loss  of  generality  we  have  replaced  dp  k  by  4  for 
convenience.  The  conditions  on  g(x)  are  then 


g(l-)  •  0-,  Real  (g)  <  0  for  0  <  X  <  1 »  (D4.D5) 


the  former  to  match  with  (Dl).  The  required  solution  of  (D3)  can  be  written  in 
the  parametric  form 
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g 


X3(v2-  v)  with  x12  ■  v'3(4  -  3v)-3  . 


(D6.D7) 


Hence  g  remains  real  and  negative  for  the  range  1  >  x  >  Xo  *  5~ 5/12  22/3 
(»  0.41516),  or  1  >  v  >  1/2,  with  (D4)  satisfied  as  v  +  1-,  x  ♦  1-.  As  v  ♦  l/2+, 
dx/dv  ♦  0  and  so  then  d2g/dx2  has  an  inverse  square  root  singularity  although 
g(<  0),  dg/dx  (>  0)  remain  finite.  The  singularity  in  g  is  accompanied  by  one  in 
A(x)  but  these  are  smoothed  out  in  a  thinner  layer  near  X  *  Xo  by  ®eans  of  Airy 
functions.  Consequently  g  becomes  complex  in  x  <  Xo*  However,  (D6) ,  (D7 )  can 
continue  to  hold  then,  with  v  now  complex,  and  the  constraint  (D5)  is  maintained 
for  0  <  x  <  Xo*  x  *  0+  instance  we  obtain 


g(x)  ~  3“5/4e'5lTi/8  +  0(x3/2) 


(D8) 


from  (D6),  since  v  ~  3  5/8x~3/2  eXp  C  — 5tt i / 16 )  from  (D7),  and  (D8)  satisfies  (D5). 


It  follows  that  the  decay  or  upstream  influence  for  small  positive  x  is 
purely  exponential  in  1  >  x  >  Xo  ^ut  oscillatory-exponential  in  Xo  ^  X  >  0. 

An  extra  viscous  layer  also  arises  nearer  the  surface  but  the  Airy  equation  holds 
again  there  and  accommodates  the  exponential  growth  required  by  (D8)  at  the  upper 
extremes  of  the  viscous  layer. 
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Figure  1  Nonuniqueness  of  Neer-Weke  Solutions 
(b)  Velocity  Function 


83-3-61-3 
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I**-”  TRIPLE  DECK  REGION 


Figure  1  Nonuniqueness  of  Neer-Wske  Solutions 
(c)  Implied  Flow  Pettems 
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Figure  1  Nonuniqusnsss  of  Ntar-Wake  Solutions 

(d)  Very-Near-Wake  Displacsmsnt  Distribution 
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G'(o)  =. 


Figure  2  More  General  Near  Wake 

(a)  Velocity  Function  Distribution 
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G'(o)=  +1 
G*(o)=  -1 


Figure  6  Eddy  Closure  on  a  Solid  Surface 
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J 


q.  REAL  (q) 
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Figure  B3  Limit  Form  Exponent 
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